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ABSTRACT: We compute the drag force experienced by a heavy quark that moves 
through plasma in a gauge theory whose dual description involves arbitrary metric and dila- 
ton fields. As a concrete application, we consider the cascading gauge theory at temperatures 
high above the deconfining scale, where we obtain a drag force with a non-trivial velocity 
dependence. We compare our results with the jet-quenching parameter for the same theory, 
and find qualitative agreement between the two approaches. Conversely, we calculate the jet- 
quenching parameter for A/" = 4 super- Yang-Mills with an R-charge density (or equivalently, 
a chemical potential), and compare our result with the corresponding drag force. 
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1 Introduction and Summary 



Recently there has been a surge of interest in the possibihty of employing the gauge/gravity 
duality [H [2] to determine the rate of energy loss in finite-temperature strongly-coupled 
gauge theories. The ultimate aim of this program would be to make contact with current 
[3] and future ^ experimental studies of quark-gluon plasma (QGP), but at this point the 
gravity dual of QCD is not yet available, so one must still be cautious when attempting to 
draw inferences in this direction. 

In the AdS/CFT context, the issue of energy loss has been approached from three different 
perspectives. The authors of [5] proposed a model-independent, non-perturbative definition 
of the jet-quenching parameter (which in the QGP case codifies the suppression of hadrons 
with high transverse momenta) in terms of a light-like Wilson loop, which they then com- 
puted in A/" = 4 super- Yang-Mills (SYM) using the recipe of [6]. Their computation was 
generalized in [7j to the cascading gauge theory [HI [9] at finite temperature [10] , and in [TT] to 
certain marginal deformations of the JV = A theory. Related work may be found in [121 [13] • 

A second approach was pursued in [Ml |T5] , where the drag force experienced by a heavy 
quark that moves through A/" = 4 SYM plasma was determined by considering a string in 
the dual AdS-Schwarzschild geometry. This calculation was extended to all asymptotically 
AdS geometries in [16], including the case dual to A/" = 4 SYM with a non-zero chemical 
potential, a case that was studied simultaneously in |17] from a different perspective: whereas 
[T7] worked directly with the ten-dimensional spinning D3-brane background, [16] employed 
instead the five- dimensional charged black hole solution of A/" = 8 gauged supergravity [H] 
obtained upon Kaluza-Klein reduction on the [I9|, thereby arriving at different results. 
A more detailed picture of the energy- loss process was painted in [20j, which studied the 
wake left by the quark as it ploughs through the plasma, using the methods of [21], |22]- The 
connection between the rate of energy loss found in [TH |15] and magnetic confinement was 
explored very recently in |23j . 

A third approach [21] extracted the diffusion coefficient for a heavy quark in the A/" = 4 
plasma from an analysis of small fluctuations of a Wilson line that follows the Schwinger- 
Keldysh contour. On the AdS side of the duality, this involved a study of fluctuations that 
propagate along a string; similar calculations were carried out in [14J. 

It is important to explore the relation between these three approaches. As a step in this 
direction, in this paper we compare the drag forces and jet-quenching parameters for two 
different gauge theories. We begin in Section [2] by generalizing the drag force calculation of 
|14[ [T5[ [16] to backgrounds with arbitrary metric and dilaton fields, finding the general result 
f l2.12p . In Section [3] we then specialize to the cascading gauge theory, where the resulting 
drag force (13.41) is found to display a highly non-trivial velocity dependence. We compare 
this force with the jet-quenching parametei0 (13.91) determined in [7], finding agreement in 

^To be more precise, we find that the directly comparable quantities are the jet-quenching parameter q 
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functional form, and numerical agreement (up to an overall constant) in the region of large 
velocities. In Section H] we proceed in the opposite direction, computing the jet-quenching 
parameter for the Af = 4 plasma with an R-charge density, and comparing the result with 
the drag force determined in ^7\. This comparison is interesting because in this case the 
quantities to be compared are not just numbers but functions of the charge density J. For 
small values of J, our result (14.151) is again in qualitative agreement with the drag force 
(14. 16p obtained in [T7], but for arbitrary charges, the results (14.121) and (I4.19P disagree. 

The general lesson appears to be that the parameter q defined by [5] and the dissipative 
force extracted from the procedure pioneered in [HI |T5] represent closely related but not 
identical measures of the rate of energy dissipation in a given non-Abelian plasma. Our 
results underline the interest in exploring the connection between the various approaches to 
energy loss from a more general viewpoint, by attempting to extrapolate from one to the 
other directly at the level of the corresponding AdS / CFT calculationsI§ 

2 Drag Force in Gauge Theories with Holographic 
Duals 

Consider a background dilaton field 0(x) and stationary Einstein frame metric 



that holographically encode the dynamics of a strongly-coupled gauge theory. Since we are 
interested in studying this gauge theory at a finite temperature T, we assume the geometry 
(14. ip includes a black hole [28j. In this setup one can introduce an external quark in the 
gauge theory by considering a string that has a single endpoint at the boundary and extends 
all the way down to the horizon [6] (the gauge theory is therefore non-confining at the given 
temperature). 

The relevant string dynamics are captured by the Nambu-Goto action 



with QafB = G^ydaX^^dpX" . Letting z denote the radial coordinate of the black hole geometry 
and t, (i = 1, 2, 3) label the directions along the boundary at spatial infinity, we make the 
static gauge choice a = 2;, r = t, and, following [TU [15], focus attention on the configuration 



and the ratio /x/T, where ^ is the friction coefficient that determines the drag force through dp/ dt = — /xp. 

^After the first version of this paper had appeared on the arXiv, the relation between the drag force 
|14|, 115] and jet-quenching [5 approaches was discussed in the context of a study of mesons moving through 
the plasma [25l [26] . Related work may be found in [27] . 



ds\ = Gi^i_„{x)dx^dx 



(2.1) 




(2.2) 



X\z,t) =vt + p{z) , X^ = = X\ 



(2.3) 
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For the appropriate sign of p', this describes the string traihng behind its boundary endpoint 
as it moves at constant speed v in the direction, a configuration dual to the external quark 
traversing the plasma. 

Using fl2.3l) in fl4.5p . we find the Lagrangian 

C = e^/^^/^d^ = e^/V-G.^Gii - G^^Cv^ - CGup'^ , (2.4) 

which results in an equation of motion for p implying that 

7rx = ^ = e'f/'^^p' (2.5) 
dp' V-^ 

is a constant. Inverting this relation we obtain 



^ GxxGtti^^GxxGtt + 7r|^) 



Just like in the TV = 4 case analyzed in [T5] , for v"^ > the numerator in this expression 
changes sign at a radius z = defined by 

iGu + G,y)U=,^ = , (2.7) 

and so the string will not be able to extend all the way down to the horizon at z = zh unless 
the denominator also changes sign at Zy. This condition fixes 



TT 



2 - -e^GxAU=.. = v'e^GlU=z^ . (2.8) 



X ^ '-^xx^tt\z=ZT] ^ ^ ^xx\z=Zv 



With (12.61) we can then compute the current density for momentum along x^, 

P: = -^.e'^^'G^^g'^d^X'' (2.9) 



1 ^/2^xxGu / 



and use it to compute the drag force experienced by the string/ quark, 

I — T^DZ s/2^xxGtt , 

— = v/=dii^P. = -^f , (2.11) 

which after some algebra is easily seen to simplify to 

^ = -I^ = -^e^/'G I (2 12) 

This generalizes the result obtained in [16] to backgrounds with an arbitrary dilaton profile. 
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3 Drag Force in the Cascading Plasma 



Let us now apply the results of the previous section to an interesting concrete example: the 
Klebanov-Strassler cascading gauge theory P [H] , whose dual geometry at temperatures high 
above the deconfining transition was constructed in [10] (see also [33|, [32]) and is given by 



[1 - z)dt^ + dX^) + 



32 



*^-2p2(r,-50. 



K+el + el+el) 



V 



2z + [-2z + {z-2) log(l - z)] log z + {z- 2)Li2(z) 

AOK^z 

o 

[log2;log(l - z) + U2{z)] , 



16K^z 



1. 



— -— + -Li2(l-^) 



24 4 



where the radial coordinate z runs from the horizon at 2; = 1 to the boundary at z 
Using (13.11) in (12. 8p one finds that z^ = 1 — v"^ and 




,2P2,, ,^/2 



or, more explicitly, 

TTx = - 
f{v) = 7 




exp 



p2 

(l + i;2) 



-Llot' — , 

24 4 4(1 -t;2) 



The drag force is then given by (12.121) as 

1 



dp 

dt KT 




27r/2 V 



(log(l-t;2)logt;2 + Li2(l-t;')) 



p2 / p4 



(3.1) 



0. 



(3.2) 



(3.3) 



(3.4) 



where we have kept only the first two terms in an expansion in powers oi P^/K^, <^ 1, 
because the solution (13.11) itself is only valid to this order. The velocity-dependence seen in 
the first term is just the p/m factor present already in A/" = 4 SYM [HI |15]. The second 
term has an additional non-trivial dependence codified in the function f{v) defined in (13. 3p . 
which approaches the value — 7r^/24 ~ —0.411 for f — 0, and has a logarithmic divergence 
in the ultrarelativistic region f — 1. As seen in Fig. 1, away from this region f{v) is nearly 
constant. 
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Figure 1: The function defined in (j3.3p . which codifies the dependence of the drag force (j3.4p on 
the velocity beyond the expected relativistic dependence v/y/l — v"^ found already in = 4 SYM. 



As in [7], it is instructive to compare the result (13. 4p for the cascading plasma against the 
drag force in A/" = 4 SYM [H I 



dp 

dtM=A 



by considering the ratio 



dp 

dtKT 
dp 

dtAf=4: 



la 1 



1 



(3.5) 



(3.6) 



Employing the relation 8a = T^Kl/A [lOj, this reads 

1 



dp 

dtKT 
dp 

dtJV=4: 



p2 / p4 



(3.7) 



Upon synchronizing the rank of the A/" = 4 theory with the effective rank of the cascading 
gauge theory at the given temperature by setting K^, = = VAk^tiN = 2^7T^gslgN, 

and using = 4:Trgs, we are left with the final result 



m p2 / p4 \ 



dp 

dtA/'=4 



(3.^ 



This ratio has the same qualitative form as the one computed for the corresponding jet- 
quenching parameter calculated in [7], 



Qcascade , P , f ^ \ 



(3.9) 



with X — —1-388. In both cases the temperature-dependence arises only from the dependence 
of the effective rank K^: on T, the precise form of which is i^'*/2P^ ^ log(T/A) [TO]. In 
addition, f{v) is negative for any value of v (see Fig. [1]), which implies that, just like the 
jet-quenching parameter, the drag force increases with increasing temperature. It also seems 
worth pointing out that precise numerical agreement between f{v) and x is achieved at a 
rather large value of the velocity, v ~ .994, which appears to be related to the fact that the 
jet-quenching calculation focuses on ultra-relativistic quarks. 

Before closing this section we should note that it has recently been argued [31] that (13.11) 
cannot be trusted all the way down to the boundary at 2; = 0, because the perturbative 
expansion in P"^ jK^ <^ 1 through which it was derived [in] breaks down at a critical radius 
Zc = e~^'^*/^ , which is non-perturbatively small but non-zero. Since the drag force (13.41) 
is determined by the value of the metric and dilaton at z^, = 1 — f^, this does not affect 
our result as long as Zy 2> Zc, which means that we can consider ultra-relativistic velocities 
except for a region parametrically close to f = 1, defined by the condition ln7 > K^^/P"^. As 
an example, validity of (13.41) at the velocity v = 0.994 considered in the preceding paragraph 
merely requires that P^/K^, ^ 0.45. The problem of obtaining a solution valid all the way 
down to 2 = has been addressed numerically in [Mj. In any case, our main goal in this 
section has been to obtain a drag force comparable to the jet-quenching parameter (13. 9p . 
which was derived in [7] using the background (13. ip . 

4 Jet-Quenching Parameter in a Charged J\f = ^ Plasma 

The near-horizon metric for rotating D3-branes at finite temperature and with one angular 
momentum turned on is [SQl [291 EI] 



1 \{^-h) 



{{dx^Y + (dx'Y) - (1 + h)dx'^dx~ + dxl + dxl + 



(4.1) 




2 




sin^ Odtdcj) + r'^{Md'^ + A sin^ Odcj)^ + cos' 



?ddnl) , 
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where 



™± _ ^ (^0 _|_ 

v2 



H 



A = 1 + 



A = 1 + 



/i 



r^A ' 

i?^ = 47riV(7,/^ (4.2) 

This background is dual to TV = 4 SYM theory with an R-charge density, or equivalently, a 
chemical potential. The geometry has an event horizon at the positive root of h{rH) = 0, 



(4.3) 



and an associated Hawking temperature, angular momentum density, and angular velocity 
at the horizon 

which translate respectively into the temperature, R-charge density and R-charge chemical 
potential in the gauge theory. 

In this section we will calculate the jet-quenching parameter qj following [3], and compare 
with the drag force result of [17J. For this we must consider a string whose endpoints lie 
at r ^ oo in the spacetime (14.11) and trace a rectangular light-like Wilson loop of length 
L along = y and along x^ [5j . Making the static gauge choice a = y,T = x^, the 
relevant configuration is r{y,x~) = r{y), with all other embedding fields constant, and with 
boundary conditions r(±L/2) = oo. The Nambu-Goto action reduces to 




S = ^KZ7r dyJ[- + ^ Ul-h). (4.5) 



Notice that, just like in [5], [TIE], we are working in Lorentzian signature and have omitted 
a factor of i in (14.51) . 
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Regarding y as 'time', the fact that the Lagrangian is time-independent imphes that the 
'energy' is conserved, a statement that is easily seen to lead to 



with 7 > R^/tq an integration constant. It follows from this equation that the minimum 
value of r(y), which by symmetry must lie at ?/ = 0, occurs at the radius where h = 0, i.e., 
at the horizon rn- 

Integrating (14.61) we find a relation between 7 and the width L of the Wilson loop, 

L 1 I"" ^ [H R^I 

dr\ T- = ' (4-7) 



(4.^ 



2 ^M/R^ - 1) Jr„ V /i thVM/R' 
where we have defined 

dp r dC 



Also, using (14. 6 p in (14.51) we are left with a trivial integral that yields 

According to the recipe of P, to compute the Wilson loop we must subtract from S the 
self-interaction of the isolated quark and the isolated antiquark, which in the AdS side 
corresponds to the Nambu-Goto action evaluated for strings that extend from (to) r ^ 00 
to (from) r = th a,t fixed 

s„=^r*,/iHi=^i!i^. (4.10) 

" 2„a' 7,. V ft 2!raV„ ^ ' 

It is clear from (14. 7p that small L corresponds to large 7. In this regime, the leading term in 
(14. 7p implies ^7 oc 1/L, which when substituted in (14. 9 p gives an L-independent result that 
is precisely cancelled by (I4.10p . The quantity we are after comes from the next-to-leading 
term in (14. 7p for large 7, which yields 

S, = S-So= . (4.11) 

Using the definition of [5], the jet-quenching parameter then follows a^ 



We employ here the final normalization of q given in v3 of [5 ; the original definition was a factor of \/2 
smaller. 



8 



The final step would be to express the result (14.121) in terms of gauge theory quantities. 
This can be done analytically in the / <^ tq regime, where up to OII'^/tq) corrections the 
relations (14. 3 p and (14. 4p imply 



4r, 



— ) = ^0 ( 1 



4J^ 



-kR^T 1 + 



4J2 



(4.13) 



and the integral (14. 8 p can be seen to give 
, v^r(5/4) 



r(3/4) +¥l(^(-^)-^^^-^))-^-4.r 



(4.14) 



where E and K respectively denote complete elliptic integrals of the second and first kind. 
The numerical value of the coefficients defined above is a ~ 1.311,6 ^ 0.7120. Putting this 
all together, and remembering that the 't Hooft coupling is given by A = Qym^ ~ R^/(^' ■, 
we are finally left with 



7r 



(4.15) 



As a (rather mild) check, note that the leading term, gj=o, reproduces the result of [5j. The 
next-to-leading term is a new result. 

Again, it is natural to compare (14.151) against the drag force calculated in [17] (see also 



dpi 
dt 



2 m 



1 + 



+ 0{J^/T 



12\ 



(4.16) 



The comparison is especially interesting because the quantities to be compared are now 
functions of the additional parameter J. The leading term in (14.161) is of course the result 
computed in [HI [15] at zero chemical potential. As in the cascading gauge theory case 
analyzed in the previous section, we find that the subleading terms in the jet-quenching 
ratio 
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l + ;^(l + V2a) + 0(^) 



and the drag force ratio. 



{dpi/dt)j 
{dpi/dt)o 



1 + 



8J2 



0{J^/T 



12\ 



(4.17) 



(4.18) 



have the same qualitative form. The numerical coefficients are also in rough agreement: they 
are both positive, implying that both the drag force and the jet-quenching parameter increase 
with increasing charge density, and they are of the same order of magnitude. 
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It is important to note, however, that this agreement cannot persist at arbitrarily high 
order in the expansion in powers of //ro, because the full result f l4.12p for the jet-quenching 
parameter evidently has a different functional dependence on ro and I than the full drag 
force result |T7j 



codify somewhat different information on the process of energy loss in a plasma. 

After the first version of this paper had been posted on the arXiv, three other calcula- 
tions of g in a charged A/" = 4 SYM plasma appeared [35l [36l E?]. The first of these is 
not directly comparable to ours, because the authors of [35] employed a five- dimensional su- 
pergravity perspective instead of the ten-dimensional string theory viewpoint adopted here 
(for a discussion on the relation between these two approaches, see [17j). Our full result 
f l4.12p for the singly-charged plasma agrees with the one obtained by the authors of [36] 
(who also determined q for two equal non-zero charges) and [37] (who addressed the general 
three-charge case). This result was plotted in [36] and shown to be non-monotonic beyond 
the small-charge region explored in (14.151) . which motivated the authors of [17J to generate 
a comparable plot of their drag force result ( 14.19^ . As expected from the discussion above, 
the two graphs are different but qualitatively similar, so again the general lesson seems to be 
that the parameter q defined by [5] and the dissipative force extracted from the procedure 
pioneered in [El |15] represent closely related but not identical measures of energy dissipa- 
tion in a given non-Abelian plasma. Evidently, more work will be required to completely 
elucidate the relation between these two approaches. 
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(4.19) 



This appears to imply that, in the general case, the drag force and jet-quenching parameter 
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